Spherical Vesicles Distorted by a Grafted Latex Bead: An Exact Solution 
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We present an exact solution to the problem of the global shape description of a spherical vesicle 
distorted by a grafted latex bead. This solution is derived by treating the nonlinearity in bending 
elasticity through the (topological) Bogomornyi decomposition technique and elastic compatibility. 
We recover the "hat-model" approximation in the limit of a small latex bead and find that the 
region antipodal to the grafted latex bead flattens. We also derive the appropriate shape equation 
using the variational principle and relevant constraints. 
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I. INTRODUCTION 

Encapsulation, binding and adsorption of particles 
onto a membrane play an important role in biological 
processes. A prominent issue in the context of vesicles is 
the understanding of long-range phenomena connected to 
the inclusion of mesoscopic particles in their membrane 
The difference in scale between the vesicle's 
membrane thickness and the long-range limit {e.g., vesi- 
cle size) allows the vesicle to be described as an embedded 
surface. The description of a vesicle's membrane behav- 
ior in the long-range limit is currently founded on the 
concept of bending elasticity Recently the Bo- 

gomol'nyi decomposition technique p^ . [Tl| has emerged 
as a promising theoretical framework to study long-range, 
nonlinear elastic phenomena in soft-condcnscd matter 

Exact solutions to nonlinear models offer novel insights 
into physical systems not at first unveiled by approximate 
approaches. Therefore our goal here is to derive and dis- 
cuss an exact solution, which is motivated by the ex- 
perimental observation of a spherical vesicle distorted by 
a grafted latex bead within the Bogomol'nyi frame- 
work combined with elastic compatibility and a global 
constraint. To this end, in the next section we intro- 
duce the bending Hamiltonian, describe the fundamen- 
tal theorem of surface theory, introduce the Bogomol'nyi 
decomposition technique, and impose a conformally in- 
variant global constraint meant to mimic the total mean 
curvature, area, and/or volume constraints. In Sec. IIIII 
we invoke the variational principle to derive the relevant 
shape equation. In Sec. IIVI we derive an exact solution 
for the deformation of the vesicle in the presence of a 
grafted latex bead at its north pole. This solution is ob- 



tained by noticing that the polar angle of the outward 
surface normal is a kink soliton solution of the double 
sine-Gordon equation. We also discuss the resultant ge- 
ometric frustration. We summarize our main findings 
in Sec. |Vl Finally, the details of the metric and shape 
tensors in isometric azimuthal coordinates are given in 
Appendix |^ whereas the shape equation is discussed in 
Appendix |B] 



II. MODEL 

A. Preliminaries 

In order to take advantage of the Bogomol'nyi tech- 
nique, we should describe the vesicle shape within a co- 
variant field theory [3, [IE]. According to the funda- 
mental theorem of surface theory [l^, [13] , any embedded 
surface is perfectly represented by a pair of symmetric 
second-rank tensors coupled to each other by integrabil- 
ity conditions: a prescribed metric tensor gy- coupled to 
a prescribed shape tensor hij via elastic compatibility 
conditions. Within this framework, the bending Hamil- 
tonian suggested by Canham Q reads: 



Tib [5] = ifc/dS by b*- 



(1) 
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which depends on the vesicle shape S through the pre- 
scribed pair (gij , hij ) and on the phcnomcnological bend- 
ing rigidity k. We have adopted the Einstein summation 
convention [l8| . and used customary notation: the inte- 
gral runs over the surface manifold S with surface ele- 
ment dS = da;^-\/|gj, where \g\ represents the determi- 
nant det(gij) and x the set of arbitrary intrinsic coor- 
dinates. Since the roundness of the grafted latex bead 
imposes axisymmetric deformations, we focus on vesi- 
cles of revolution. In a recent work [l^, we have shown 
that, for surfaces of revolution, the isometric azimuthal 
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coordinates (u, ip) are a relevant choice, which not only 
drastically simplifies formulas related to the fundamen- 
tal theorem of surface theory but also elucidates the ap- 
plication of the Bogomol'nyi technique to the bending 
Hamiltonian ([T]). 



B. Fundamental Theorem of Surface Theory 

Working with abstract manifolds described by a 
parametrized metric tensor gy is common in general rela- 
tivity [l^: when the abstract manifold is embedded, the 
metric tensor may be coupled with a shape tensor by- . 
For embedded bidimensional abstract manifolds, the cou- 
pling occurs through integrability conditions (the Gauss- 
Codazzi equations) and a system of differential equations 
to integrate (the Gauss- Weingarten equations under ap- 
propriate additional conditions): this is the substance 
on which the statement and the proof of the fundamen- 
tal theorem of surface theory are based [H, Fol- 
lowing the fundamental theorem of surface theory allows 
us to claim that any pair of diagonal second-rank ten- 
sors (gjj,bij), which in isometric azimuthal coordinates 
(u, iy9), takes the form 



„2a(«) 



s''(")9„17(u) and 



(2) 

-c'"("'sinf7(u), (3) 



and obeys the clastic compatibility condition 
du<j{u) — cosfi(w), 



(4) 



where the local Weyl gauge field a and the polar angle 
f2 of the outward surface normal are sufficiently differen- 
tiable functions of u, corresponds to a unique axisymmet- 
ric surface S — modulo its position in space. Detailed in 
Appendix \^ the demonstration consists in showing that 
the associated Gauss- Weingarten equations under appro- 
priate additional conditions determine a unique surface 
of revolution modulo an arbitrary rigid motion. In simple 
terms, straightforward successive integrations [computa- 
tional chain (|A32p - (|A50p ] reduce the set of differential 
equations to integrate [ (jA26|) - (|A28|) ] to the pair of equa- 
tions: 



r(u) =e'^("\ 
duz(u) = -e""*") sinfJ(u), 



(5) 
(6) 



where r{u) and z{u) are, respectively, the radius and 
the height of the axisymmetric surface in cylindrical 
parametrization. Formula ([5]) is a reparametrization that 
gives a posteriori a second interpretation to a. In short, 
for abstract surfaces of revolution (H))"®, the integrabil- 
ity conditions reduce to the single condition ([4]) while the 
system of differential equations to be integrated simpli- 
fies to the single equation Therefrom, the bending 
Hamiltonian ([1]) can be written as 



Tib [S] =TTk du 
Jv 



sm 



(7) 




FIG. 1: Cross-sectional profiles of a spherical vesicle of rev- 
olution distorted by a round latex bead grafted at its north 
pole with respect to different encapsulated radii p. The bold 
arcs indicate the polar cap imposed by the bead. Dimension- 
less latex bead radius p — 0.15; relative encapsulated radii p/p 
from outside to inside: 0, |, |, ^, |, |, |, and 1. 



where the integral runs along the profile V of the axisym- 
metric surface S. 



C. Bogomol'nyi Decomposition 

The Bogomol'nyi technique allows us to resolve the 
Hamiltonian ([1]) into a perfect square Hamiltonian and 
a topological bound [ij, for spherical surfaces of 

revolution So, the decomposition reads [l^ 

Tib [So] ^TTk [du [dun{u) - sin 17(u)] V 47rfc, (8) 
Jv 

which readily saturates the bound when the polar angle 
n satisfies the first-order nonlinear differential equation 



9„51(w) = sinri(u). 



(9) 



The centered solution of © is the aperiodic sine-Gordon 
kink: 



il{u) ~ 2 arctane". 



(10) 



As expected [IJ, [l5[ , the axisymmetric surface obtained 
by integrating the related Gauss- Weingarten equations 
under appropriate conditions f Appendix E]) is the round 
sphere, in cylindrical parametrization. 



r{u) — sechu, z{u) 



tanhu, 



(11) 
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the isometric coordinate u running from —oo to +00. 
This remarkable result has encouraged us [l5| to envision 
deformed spherical vesicles of revolution as frustrated or 
unsaturated sine-Gordon kinks lla. |2(]|. I21I. 



D. System of Investigation 

Now let us assume a "bare" spherical vesicle of revolu- 
tion distorted by a round latex bead chemically grafted 
at its north pole [1]. By bare, we mean that only the 
bending energy ^ is considered. Furthermore, to mimic 
the total mean curvature, area and/or volume constraints 
, we impose the following covariantlikc global con- 
straint: 



Constant. 



(12) 



As the metric determinant \g\ is not a (covariant) scalar, 
the global constraint ([T2|) is not covariant: in other 
words, because it depends on the choice of coordinates, 
the global constraint is apparently unphysical. On 
the other hand, the particular choice of the exponent ap- 
plied to the metric determinant \g\ renders the global con- 
straint (jl2p locally invariant under conformal transforma- 
tions: hence, contrary to the (physical) area and/or vol- 
ume constraints, the (unphysical) global constraint (|12[) 
is conformally invariant just like the bending Hamilto- 
nian ([T]). Also the intuitive motivation behind the choice 
of p2p as a global constraint appears naturally in iso- 
metric azimuthal coordinates (u, tp) as then we have 



27r du 
Jv 



Constant. 



(13) 



Thus, ultimately our system is governed by the Hamil- 
tonian ([7]) [or ([8])] which must be minimized subject to 
the clastic compatibility condition ([U, to the geometrical 
confinement ([H]), and to the global constraint (H 



III. SPHERICAL SHAPE EQUATION 

A. Variational Principle 

Applying the method of Lagrange multipliers for differ- 
ential equation constraints and integral constraints (23j 
gives the unconstrained functional 

1 [S] ='Kk I du \{du^{u)Y + sin^ n{u) 
Jv L J 

-I- 27rfc / du /i(u) [cosil(u) — 9u(T(7i)] 



2'nk j du v{u ) [(9„z(u) -He''(")sinf}(u)] 
Jv 

+ 2'Kk a I du, 
Jv 



(14) 



where /i and v arc local Lagrange multipliers while a is a 
global Lagrange multiplier. Whereas the Euler-Lagrange 
equations derived from the functional pi)) by varying 
with respect to /i and v reproduce, respectively, the elas- 
tic compatibility condition ([4]) and the geometrical con- 
finement ^ as required, the ones with respect to cr, Vl, 
and z arc, respectively, 



^u^i{u) = -v{u) e'^(") sinf}(u), 
duu^{u) ~ sin ri(u) cos r2(u) 

- ^{u) s\x\Q,{u) + v{u) e''^") cosfi(u), 
duv{u) = 0. (15c) 



(15a) 
(15b) 



Furthermore, since the integrand of the functional (|14p 
does not depend explicitly on the variable of integration 
u, the Beltrami identity |23| may be computed; it yields 



\{dun{u)y 



isin^ ^l{u) + /i(u) cosf2(ii) 
+ v{u) e'"(") sin Q.{u) + a ^ J, 



(16) 



with J a constant of integration. On the other hand, at 
the boundary dV, arbitrary small variations of a, 57, z, 
and u — denoted by da, SO., Sz, and 6u, respectively, — 
must hold: 



du^{u) 5fl{u) 
v{u) Sz(u) 
J Su 





= 0, 


dV 




-0, 


dV 




= 0, 


dV 




= 0. 


dV 





(17a) 
(17b) 
(17c) 
(17d) 



Clearly the conditions (|17ap and (|17b[) are always sat- 
isfied at both poles: at the south pole the radius must 
vanish (5a = 0) and the surface normal must stay parallel 
to the axis of revolution ((5ri = 0), whereas at the north 
pole the grafted latex bead dictates both a fixed radius 
(5a = 0) and a fixed polar angle (5il, = 0) along the in- 
terface parallel Let us now ignore for a while the 
condition (|17cp and focus on the last condition (jl7d|) . In 
general, because there is no restriction on the meridian 
arc-length between the podal boundaries 25| — namely 



because there is no restriction on the isometric coordinate 
u — , the variable u must be permitted to vary freely at 
the endpoints, therefore the Beltrami constant of integra- 
tion J must be set to zero in order to fulfill the condition 

(fTTdli . 

However, since the imposed global constraint (fT^ . ac- 
cording to Eq. (Uni), fixes the length of the interval along 
which the variable u must vary, the variable u is (inten- 
tionally) not allowed to vary freely at the endpoints here. 
Thus, here, the Beltrami constant of integration J is (ar- 
tificially) an arbitrary constant. Nonetheless, we have 
to remember that the imposed global constraint ([T2)) [or 
P3p ] is meant to mimic true global constraints. Besides, 
it is noteworthy that the Beltrami constant of integra- 
tion J can be absorbed by the artificial global Lagrange 
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multiplier a associated to this global constraint. Below, 
with respect to our approach so far, the Beltrami con- 
stant of integration J will be set to zero in order to meet 
the general statement associated to the condition (jl7dp 
and in such a way that any "artifact" introduced by the 
imposed constraint (fT2|) will only be carried out by the 
global Lagrange multiplier a. Meanwhile, as an immedi- 
ate consequence of Eq. (jlScp . the discerning reader has 
noticed that the local Lagrange multiplier ^{u) is a con- 
stant: 

= 1^0- (18) 

Promptly, having still in mind the geometrical confine- 
ment ([6|), the same reader has furthermore resolved 
Eq. (|15ap — formally at least: 

H{u) = vq z{u) + ^iQ, (19) 

with /io a constant of integration. Therefore, the Bel- 
trami identity pB]) writes 

- ^{dui^{u)Y + isin^ n{u) + fio cosfi(w) 

-I- 1^0 [z{u) cos il{u) + r{u) sin il(u)] + a = 0, (20) 

where the reparametrization formula ([5]) has been used. 

B. Generic Shape Equation 

As a matter of fact, for generic closed surfaces of rev- 
olution, the conditions (|17ap and (jl7bp are in general 
verified because the local Weyl gauge field a and the 
polar angle U do not vary freely at the endpoints. As 
concerns spherical surfaces of revolution, their arbitrary 
small variations at each pole have to vanish, either as 
required by continuity and symmetry in the absence of 
any contact, or as dictated along the junction interface 
in the presence of a contact. Whereas for free toroidal 
surfaces of revolution, because the boundaries coincide, 
their arbitrary small variations on each side of the com- 
mon boundary have to coincide. By now, the discern- 
ing reader has readily checked that the differentiation of 
the Beltrami identity ((20)) followed by the substitution 
of the formulas (g]), ©, dH]), and actually leads 
to Eq. (|15b|) as expected. In other words, the Beltrami 
identity (j20p is the generic shape equation of our system. 
As demonstrated in Appendix |Bl the generic shape equa- 
tion (j20p becomes the customary shape equation found 
in the literature when the "artifact" disappears. How- 
ever, the generic shape equation (j20p may be specialized 
to spherical surfaces of revolution as follows. 

C. Topological Specialization 

As observed before, the bending Hamiltonian ([T|) can 
be specialized with respect to the underlying topology 
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s 

FIG. 2: The Gaussian curvature K versus the length s along 
meridians from the south pole for different encapsulated radii 
p of the northern grafted latex bead: the encapsulated radius 
p increases in the direction of the arrows. Dimensionless la- 
tex bead radius p = 0.15; relative encapsulated radii p/p with 
respect to the arrows: 0, j, |, i, |, |, |, and 1. The inset 
shows the geometry at the north pole: the dashed circle pro- 
files the latex bead, and the bold arc the polar cap imposed 
by it; the arrows indicate the nomenclature: the latex bead 
radius p and the encapsulated radius p. 

through the Bogomol'nyi decomposition technique 
[Tsl : it appears that the generic shape equation (|20p can 
be specialized with respect to the underlying topology as 
well. Indeed, the vanishing of the constant of integration 
vq depends on the arbitrariness of the small variations of 
the height z along the axis of rotation at the boundaries 
dV, that is to say on the underlying topology. While 
for spherical surfaces of revolution the height z must be 
free to move along the axis of revolution at each pole, for 
free toroidal surfaces of revolution the small variations of 
the height z must coincide along the common boundary: 
hence in order to meet the boundary condition (|17cp . i>o 
must be zero when the topology is spherical, whereas 
vo is a priori arbitrary when the topology is toroidal. 
Henceforth, when the underlying topology is spherical, 
the Beltrami identity or the generic shape equation ^U\\ 
becomes 

- i {dun{u))^ + isin^ n{u) + fio cos n{u) +a = 0. (21) 

The reader, who still has in mind the above implemen- 
tation of the variational principle, might point out that 
this topological specialization allows us to relax the geo- 
metrical confinement ([6|) as far as only spherical vesicles 
of revolution — governed by an Hamiltonian indepen- 
dent of the height z and its derivatives — are concerned. 
This approach has been implicitly used by us to derive 
the shape equation in the illustration part of a previ- 
ous Paper — namely Ref. |T^. Whatever, the specialized 
Beltrami identity (|2ip is the announced spherical shape 



5 



equation. Of course, for more realistic systems the term 
inherited from global constraints does not reduce to a 
constant: this concern will be the main focus of our fu- 
ture investigations. As usual, boundary conditions more 
specific to the involved system allow us to establish rela- 
tionships between the constants of integration. Here, the 
boundary conditions at the south pole (fi = tt, = 0) 
give 



a 



(22) 



Therefrom, the spherical Beltrami identity (|2T|) reads 



the Euler-Lagrangc equation with respect to Q is the dou- 
ble sine-Gordon equation ([M)) : for spherical surfaces of 
revolution So, the effective Hamiltonian ([25|) decomposes 
with respect to the Bogomol'nyi technique into: 



Tieff [So] = Trk du [9„f2(u) — sinri(w)]' 



+ 2TTke du cos n{u) + ink. (26) 
Jv 



- i(a„17(w))^ + ^sin^niu) +ecosn{u) +e = 0, (23) 

which is clearly the integrated form of the double sine- 
Gordon equation [2^, [23| 



duu^{u) = sm^l{u) cosr2(u) — e smn{i 



(24) 



Let us notice that the coefficient e arises from both the 
elastic compatibility (j4|) and the imposed global con- 
straint (fT2|) . while it is a constant because the spherical 
version of the functional pi]) does not depend explicitly 
on a. As shown below, grafting a round latex bead on our 
bare spherical vesicle of revolution allows us to impose a 
nonvanishing coefficient e through the prescribed bound- 
ary conditions. Subsequently, it might also allow us to 
acquire both a local and a global insight into the compe- 
tition between bending energy, elastic compatibility, and 
global constraints through a model which fortunately is 
tractable. 



IV. LOCAL VERSUS GLOBAL 



Exact Solution 



Fortunately a suitable solution of Eq. ((^51) [or (|24p ] 
generalizing the solution (|10p exists that meets the in- 
terface conditions imposed by the grafted latex bead as 
follows. First, we easily check by direct substitution that 
the solution 



n{u\e) — 2 arctan ( ^/ cosh i / u -\- sinh J j 




satisfies the second order nonlinear differential equation 
(Pi)) with 



D. Naive Analogy 

Physically, for the surface normal, the elastic compat- 
ibility (jl]) acts as a uniform external axial field of mag- 
nitude ke: indeed for the related effective elastic Hamil- 
tonian, 



-Heff [S] 



-irk j du 
'v 



2 £ COS J7(u) 
(25) 



^ 2e e 
£ = or e ~ - — 



(28) 



whereas it reduces to the solution (|T0|) when the parame- 
ter £ vanishes. Solving the relevant Gauss- Weingartcn 
equations under appropriate conditions (Appendix [X]) 
leads, after detailed calculations, to an exact solution for 
the surface of revolution, which is the relevant general- 
ization of (fTT|) : 



t(u I e) 



z(M|e) 




and which smoothly joins the concavely bound spheri- interface parallel u = u [2J] and the parameter e obey 
cal northern (bead) cap to the vesicle surface when the ~ 

^ e < 1, r{u\e) = p, r2(u|e) = — arcsin-. 
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p and p being, respectively, the latex bead radius and the 
encapsulated radius, the isometric coordinate u varying 
from u to +00 (see Fig. [1] and the inset of Fig. [2]). Basic 
analytical considerations always ensure that the system 
pop has one unique solution which can be found numer- 
ically without difficulty. The scale and the position in 
space have been chosen with respect to the equator [gjj : 
its radius is set to unity (its scale) while its plane passes 
through the origin (its position). It can be readily ver- 
ified that the surface approaches the round sphere 
as expected, when the parameter e tends to zero. In 
fact, the interface conditions (|30p essentially require that 
at the proper place (the inequalities) the normal lines of 
the vesicle surface and of the round latex bead surface 
coincide (the equalities) 0, 

When the latex bead is very tightly bound to the vesi- 
cle — biotinylated lipids can be employed to render vesi- 
cles very sticky to streptavidin coated latex beads [1] — 
these junction conditions are reasonable at the vesicle 
scale, nevertheless at the membrane scale a more detailed 
treatment may be needed since the curvature experiences 
a discontinuity along the junctioiiparallel — they are the 
ones implicitly assumed in Rcf. [j. The relative magni- 
tude e of the effective uniform external field, according to 
the equality ([28|) and the interface conditions ([30|) , is de- 
termined by the latex bead radius p and the encapsulated 
radius p: for physically relevant latex beads, the relative 
magnitude e increases as either the relative encapsulated 
radius p/ p or the latex bead radius p increases. 



B. Hat-Model Limit 

Finally, for complete encapsulation {p = p) in the 
limit £ ^ 1, the "hat-model" approximation 0, [H, 



|29j is recovered as the axisymmetric surface (|29|) ex- 
periences complete contact of order O (e^ In^ e) with 
a catenoid along the corresponding interface paral- 
lel, which then fuses itself with the (neck) equator 

/l^ln(^SV^V After the coordinate shift 




^±^), we get 



r{v\e) = 2e coshw -I- 2e^lne coshw -I- • • • , 
z(w I e) = 1 + 2e (In ie + i) + 2£ w + In^ e - 



since r{v = 0\e) = 2e + 0(e'^\n'^ e), the limit e < 1 is 
here equivalent to the limit p <^ 1, thus it applies to small 
latex beads. More detailed analysis may be required for 
arbitrary encapsulations (0 < /5 ^ p). 



Geometrical Frustration 



Next, while keeping in mind that grafting a second 
identical latex bead antipodal to the first one (i.e., at the 
south pole) cancels the coefficient e in equation ([M)) as it 
was implicitly shown in Ref . [l^, we investigate the global 
deformation experienced by our bare spherical vesicle. 
First, for the surface of revolution ((29|) . the length s{u \ e) 
along meridians [2^ from the south pole u = +00 to the 
parallel u, the Gaussian curvature K(u \ e), and the mean 
curvature H{u \ e) along the parallel u are found to take, 
respectively, the exact forms 



s(w|e) = 2 



K{u\e) 



H{u\e) 



+2$ -c Vi- ,l,Wl-f 



I (l-e2)i(l + y^)^ VT^^cosh2,/i±f u-hsinh2 1^+' 



16 



/T+?+\/I+3? ^ V i+s 



cosh J u + Vl — sinh J - 



4 ^ ' 



'4 ^ 



1^ 

V \/T+?+\/T^ / 



cosh y u + Vl — sinh y : 
I 



(31) 



(32) 



(33) 



with $ denoting the Lerch transcendental function [30|. 
Plotting the Gaussian curvature as a function of the 
meridian arclength s from the south pole reveals, as 



clearly exhibited in Fig. [2l that the surface of our bare 
vesicle distorted by a grafted round latex bead surpris- 
ingly splits into three domains: a strong curvature gra- 
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dient region in the neighborhood of the bead; a weak 
curvature gradient zone around the equator; and an un- 
expectedly flat southern region, as the curvature amaz- 
ingly tends to zero at the south pole. On the other hand, 
the presence (or absence) of an antipodal grafted latex 
bead does not seem to radically affect the profiles around 
the grafted latex bead, as is evident by comparing Fig. [1] 
with the corresponding one in Ref. The global defor- 
mation of our bare spherical vesicle can be understood as 
follows. By imposing a concavely bound spherical cap, 
the grafted round latex bead prevents the bending Hamil- 
tonian ([7]) [or ([5])] from reaching its global minimum and 
brings into play the elastic compatibility (|4]): the for- 
mer phenomenon is revealed through the Bogomol'nyi 
decomposition technique, the latter can be mimicked by 
an effective uniform external axial field interacting with 
the surface normal. 

This key result concisely expressed in the decomposed 
effective Hamiltonian (|26p leads to a clear comprehension 
of the geometrical frustration experienced by our bare 
spherical vesicle: the surface normal attempts to shape 
our bare spherical vesicle into a round sphere to saturate 
that part of the effective Hamiltonian (pS]) which is inher- 
ited from the bending Hamiltonian, whereas the effective 
uniform external axial field tries to antialign the surface 
normal to minimize that part of the effective Hamilto- 
nian (|26p which mimics a uniform external axial interac- 
tion. Ultimately the alteration due to the presence (or 
absence) of an antipodal grafted latex bead resolves the 
competition between the two parts: the highly distorted 
region in the neighborhood of the grafted bead is mainly 
governed by the shape part which tends to make it round; 
by contrast, the essentially flat region antipodal to the 
grafted bead is driven by the effective uniform external 
axial field which flattens it by forcing the local curva- 
ture to vanish while the shape part tries to maintain it 
round. The elastic compatibility condition ^ globally 
"releases" the bending Hamiltonian ([7]) [or ([5])] in the 
sense that omitting the elastic compatibility condition 
(PI is tantamount to enhancing the bending energy by 
attaching an antipodal latex bead. 



V. CONCLUSION 

In conclusion, we have demonstrated that bare spher- 
ical vesicles distorted by a grafted latex bead might be 
treated by taking into account nonlinearity in the bend- 
ing elasticity through the Bogomol'nyi decomposition 
technique and elastic compatibility. We have obtained 
and analyzed a suitable exact solution, which recovers the 
hat-model approximation 0, in the limit of small la- 
tex beads, through an effective field theory. The present 
study provides a physically motivated example of pure 
geometrical frustration induced by a mesoscopic parti- 
cle. Finally, whereas mesoscopic inclusion problems are 
typically treated by considering spherical biological vesi- 
cles as infinite planes 0, 0, [H, [2^, [3l| , we have shown 



that simple topological arguments may allow us to tackle 
such problems by emphasizing the underlying spherical 
topology. We have also demonstrated how our approach 
can recover the relevant shape equation. Nevertheless, 
in order to overcome one drawback of our exact solution 
— viz., the use of an unphysical global constraint which 
is meant to mimic the usual physical global constraints 
— , our suggested approach may be further augmented by 
accounting for the physical global constraints — namely, 
the total mean curvature, area, and volume constraints 
[13 — Slid, eventually, by connecting it to other geo- 
metrical approaches d, [33, [H, [13, . Also, by applying 
it to vesicle systems exhibiting pertinent global distor- 
tions — such as toroidal vesicles [s^, vesicle systems 
exhibiting domains [s^l; vesicles distorted by adsorbed 
ob j ects M , c y lindrical and other inclusions in vesicles 
[ssl. [39I. [40l. l4ll. [4^. l43j . and so forth — . our suggested ap- 
proach might lead to more valuable solutions, but they 
are difficult to find analytically. 

Partial support to J. B. from EC under Contract No. 
HPRN-CR-1999-00163 (LOCNET network) and from the 
National Center for Theoretical Sciences, Taiwan, are ac- 
knowledged. The work at Los Alamos was supported by 
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APPENDIX A: SURFACES OF REVOLUTION IN 
ISOMETRIC AZIMUTHAL COORDINATES 

1. Metric and Shape Tensors 

Let us represent any surface of revolution in azimuthal 
coordinates (u, if) by 

X(u, (/s) = ^r{v) cos ip,r{v) sin ip, z{v)^ (A-l) 

where the radius r{v) [0 ^ ''(w)] and the height z{v) are 
both sufficiently differentiable functions of the coordinate 
V. Then the metric tensor gy has the following compo- 
nents [ll,[i3|: 

gl,„ a„X(u, If) ■ 5t,X(v, (fi) 

= {d,r{v)y + {d.,z{v)y, 

gv<p = g^v = di,X{v, If) ■ d^X{v, ip) = 0, (A2b) 
g^^ = d^X{v, ip) ■ d^X{v, if) = r(u)2. (A2c) 

Before computing the shape tensor by , the outward sur- 
face unit normal vector lSI{v,ip) may be computed as 
[ll,[ll 

^ dyXjv, Lp) X d^Xjy, if) 
||3t,X(u,(/3) X 9^X(w,(/?)|| 

= , ^ (A3) 

{d,r{v)f + {d,z{v)f 

■ \—dvz{v) cos p, —dvz{v) sin ip, 9„r(v)] . 



Hence M, 



bvv = dvv^{v, If) ■ N(w, ip) 

_ dyyz{v)dyr{v) - dyyr{v)d^z{v) (A4a) 

{dyr{v)Y + {dyZ{v)y 

= by,„ = ai,^X(w, • N(w, = 0, (A4b) 



riv^dyziv) 



(a„r(«))' + (a„z(z;))' 



(A4c) 



Next we may introduce the polar angle f2(t;) associated 
with N(w, </?): 

= arctan(— 9i,2;(w), i9t,r(u)), (A5) 
which satisfies the (useful) relationships: 

dyr[v) 



cosf2(?;) 
sin f2(w) 



v/(a„r(«))' + (a„z(t;))' 



(9„r(«))' + (a„z(i;))'' 



(A6a) 
(A6b) 



^ dyyr{v)dyZ{v)-dyyZ{v)d r{v)^ 

Thus, whereas the outward surface unit normal vector 
'N{v,ip) takes the desired form 

N{v,ip) = [sin Q{v) cos tp, sin n{v) sin cos , (A7) 

the components of the shape tensor by are readily writ- 
ten as 



b™ = -^Z + {dyz{v)f dMv), (A8a) 
bt,^ = h^v = 0, (A8b) 
b^jcp = — sinri(w). (A8c) 



Now let us transform to isometric coordinates (u, (p) [la . 
by the transformation 



J r(w) 

which yields a conformally flat metric [i^l- As a matter 
of fact, the components of the metric tensor become 



0, 



(AlOa) 
(AlOb) 
(AlOc) 



while the components of the shape tensor attain the form 



bip« = 0, 



bi^ip — — 7'(u) sinil(u). 



(Alia) 
(AUb) 
(Allc) 



Finally, to emphasize the conformally flat nature of the 
metric and/or the isometric nature of the coordinates 
(m, gg), we may introduce the local Weyl gauge field a{u) 



a{u) = ln(r(u)). 



(A12) 



To summarize, for any surface of revolution the metric 
tensor gy and the shape tensor by can take, in (isomet- 
ric) azimuthal coordinates (w, <p), the form 



,2a(«) 



g'"(")a„17(u) and \ 



(A13) 
3'"(")sinl7(u), (A14) 



respectively. 



2. Elastic Compatibility Conditions 

Next, with the perspective of computing the elas- 
tic compatibility conditions (the Gauss-Codazzi-Peterson 
equations) [l3| we may compute the Christoffcl symbols 
of the second kind T^,- [l7l.l4il. 



r«.. = 

U 2 



1 „kl 



diSij + djgH - dig^j 
and then the Riemann tensor Riju p7l. |44|| . 



(A15) 



R 



ijkl 



■ ^ kn^ Ij 



-pm pn 
In kj 



(A16) 

In isometric azimuthal coordinates (u,ip), the non- 
vanishing Christoffcl symbols of the second kind V'^^j vcr- 

ify 



r\J") = -r\^iu) = r\^iu) = r^J") = 5„a(«), 

(A17) 

and the nonvanishing component of the Riemann tensor 



R — -e^""*^"^ d a(y) 
Therefore the Gauss equations [l3| 

Ryfei = bifcbj/ — biihjk, 
simplify to the equation 

duu<^{u) ~ — sinr2(u) du^iu), 
while the equations of Codazzi and Peterson [l] 
Sfcby 

yield 

du<j{u) = cosfi(u). 



(A18) 
(A19) 
(A20) 
(A21) 
(A22) 



which clearly satisfies the Gauss equation (jA20p . In 
short, the elastic compatibility conditions (|A19p - (jA2ip 
associated with the metric tensor (jA13l) coupled to the 
shape tensor (|A14p reduce to Eq. (|A22|) . 
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3. Fundamental Theorem of Surface Theory 

Next we shall show the converse of the previous re- 
sults, namely, that any pair of diagonal second-rank ten- 
sors {gij , bij ) which in (isometric) azimuthal coordinates 
{u,(p) takes the form (jA13p - (jA14p and obeys the elas- 
tic compatibility condition (|A22p . where the local Weyl 
gauge field a and the polar angle f2 of the outward sur- 
face normal are sufficiently differentiable functions of u, 
corresponds to a unique axisymmetric surface — modulo 
its position in space. In fact, this converse theorem is 
rather a simple illustration of the fundamental theorem 
of surface theory [l^, |T3] ■ The general demonstration of 
the fundamental theorem of surface theory [l6| consists 
in proving that the Gauss surface equations |l6. [TtI 



9,„X + b,„-N, 



combined with the Weingarten equations [T^ [l3| 

under the additional conditions [l^ 

= 1, • N = 0, 



9yX-N = b,,. 



(A23) 



(A24) 



(A25) 



determine a unique surface X — modulo its position in 
space — with N as its outward surface unit normal vec- 
tor. In our case, the Gauss surface equations (|A23p re- 
duce to 



■cosn{u) duX{u,(f) 

-e"(")9„17(u)N(u,(^), 
- cos fl{u) d^pX{u, Lp), 
cosn(u) duX{u, Lp) 

-e"^"^ smn{u)N{u,ip), 
and the Weingarten equations (jA24[) to 

duNiu, cp) = c"'^("'a„r!(?/) duX{u, cp), 

a<^N(w, ^) c^'^t"' smn{u)d^X{u,ip), 
while the additional conditions (jA25|) yield 

N\u,^) 
duX{u, (fi) ■ N(u, (p) 
d^X{u, ip) ■ N(w, ip) 
duX{u, p) ■ duX{u, ip) 
duX{u, ip) ■ d;pX{u, ip) 
d^X{u, ip) ■ d^X{u, ip) 
duuX{u, ip) ■ N(w, ip) 
du^Xiu, ip) ■ N(u, ip) 
9^^X(w, ip) ■ N(u, ip) 



= 1, 
= 0, 
= 0, 

- p2<T(«) 



,2<t(«) 



0, 



e'"(") sinO(u) 



(A26a) 
(A26b) 
(A26c) 



(A27a) 
(A27b) 

(A28a) 
(A28b) 
(A28c) 
(A28d) 
(A28e) 
(A28f) 
(A28g) 
(A28h) 
(A28i) 



Note that, provided that conditions (|A28a|) - (|A28cp 
are met, conditions ( |A28gP -( |A28ip are automatically 
satisfied since then they straightaway follow from 
Eqs. (|A26a[) -( jA26c[) . respectively. By definition, on the 
other hand, the outward surface unit normal vector 
N(u, ip) expresses as 



N{u,ip) 



duX{u, ip) X 9yX(u, ip) 



(A29) 



\duX{u,ip) X d^X{u,ip)\\ ' 
which; with respect to conditions (jA28dp - (|A28fj) . writes 
N(m, ip) = e-^-^f") duX{u, ip) X d^X{u, ip), (A30) 
as we find 

[duX{u, if) X d^X{u, ip)] ' = e4"(") . (A31) 



Conditions (|A28ap - (|A28cp are clearly fulfilled. In short, 
only the additional conditions (jA28dp - ()A28f| may be 
pertinent in due course. By elimination of N(m, ip), Eqs. 
(lA26cll . (|A26bll . and (|A27bll hold: 



d^^^X{u, ip) -f d^X{u, ip) = 0, 



(A32) 



X{u, ip) = cos ip K{u) + sin B(w) -f C(u). (A33) 

Then equation (jA26bp gives 

duA{u) ~ (?„cr(u) A(u) and 9uB(u) = ducr{u) B(m), 

(A34) 



A(u) = e'"("' A and B(u) e'^^"' B, 



(A35) 



where A and B are constant vectors of integration. Also, 
formula (|A33|) writes 

X{u,if) =e'^(")[cos(^A + sin</jB] +C{u). (A36) 



Moreover, eliminating N(u, ip) between Eqs. (jA26aP and 
(jA26cp . we obtain 

sinri(u) duuX{u, ip) - du^{u) d^^pX^u, ip) 

= cosf2(u) (sinri(u) + c)„0(u)) duX[u,ip). (A37) 

Substitution of formula (|A36p into Eq. (|A37j) yields 

duuC{u) = [du(j{u) + cot VL{u)duVL{u)] duC{u), (A38) 



duC{u) 



sin Vl{u) C, 



(A39) 



where C is a constant vector of integration. So, we read- 
ily have 



/u 
dw e'"("') sinf^H C-HD, 



(A40) 
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with D a constant vector of integration. Henceforth, for- 
mula (|A36p becomes 

X(it, ip) = e""*") [cos 1^ A + sin ip B] 

- dw e""*'") sin n{w) C + D. (A41) 

We must next choose the constant vectors of integration 
A, B, C, and D so as to satisfy the pertinent additional 
conditions (|A28dp - (U28^ . First, from condition jJM^ 
we obtain 



- lp2<T(«) 



(A^+B^) +cos2ip (A^-B^) - sin2^ A-B 



or 



A^ = B^ = 1 and A • B = 0. 

Therefrom, condition (|A28e[) holds: 

9„X(u, (f) ■ d^X{u, f) = 

= e2'"(") sin n (u) [sin A • C - cos B • C] , 



or 



(A42) 



(A43) 



(A44) 



(A45) 



A • C = B • C = 0. 

Then, condition (jA28dp reads 

9„X(m, ip) ■ (9„X(u, ip) = e2'^(") = e2'^(") C^, (A46) 



= 1. 



(A47) 



In other words, the three constant vectors of integration 
A, B, and C must form an orthonormal triplet of con- 
stant vectors. Henceforth, since no additional condition 
has to be fulfilled, the constant vector of integration D is 
really an arbitrary constant vector, namely an arbitrary 
translation. In order to determine the orientation of the 
orthonormal triplet, we may explicitly compute the out- 
ward surface unit normal vector N(w, ip): 

N(it, ip) = sin f7(it) cos B X C 

+ sin rj(it) sinv? C X A + cosf2(u) A x B. (A48) 



By inserting previous formula (|A48p and formula (jA4ip . 
Eq. ((A27b)) holds: 



C = A X B, 



(A49) 



thus the orthonormal triplet is direct. To stress the na- 
ture of the direct orthonormal triplet [A, B, C] and of the 
translational vector D, we may denote them by [ex, By, e^] 
and t, respectively: eventually formula (|A41[) may read 

X(u, ip) ~ c"^"^ [cos pex + sin ip Gy] 

- / dw c'"("'' sin n{w) + t, (A50) 



which is the equation of an axisymmetric surface revolv- 
ing around the axis shifted along t with radius r(u) 
and height z{u) given by 



r(w) =e'"("\ 



(A51a) 



and 



/u 
dwe'''^'"'>smn{w). (A51b) 

By the choice of the direct orthonormal triplet [gx, e^, e^] 
and of the translation t, the surface of revolution (|A50[) 
can be placed in any position in space. Finally let us 
notice that the ultimate differential equation to integrate 
(|A39p translates within this choice of notation into 



duz{u) =-e'"(") sinn{u) 



(A52) 



Curvatures 



Whereas the intrinsic curvature (or the Gaussian cur- 
vature) K characterizes the metric tensor gij through the 
Riemann tensor Rijki, 



ikjl 1 



(A53) 



the mean curvature H and the extrinsic curvature G 
characterize the shape tensor by-: the mean curvature 
H is its half trace. 



(A54) 



while the extrinsic curvature G is its determinant. 



ikDjl- 



(A55) 



The tensor e^- being the totally antisymmetric tensor: 



g\ e,j with e,j ^ [_? +J] 



(A56) 



Therefore, in isometric azimuthal coordinates (u,ip)^ the 
intrinsic curvature K holds: 



K{u) ^ -e-2-(") duu^{u), 



(A57) 



according to formula (jAlSp . On the other hand, the mean 
curvature H yields 



H{u)^~\ e-"'-''' [dun{u) + sin (w)] , (A58) 
and the extrinsic curvature G verifies 

G{u) = e-2'^(") sin n{u)dufl{u). (A59) 



By invoking the elastic compatibility condition (|A22p . it 
is easily checked that the intrinsic curvature K and the 
extrinsic curvature G are effectively the same as asserted 
by the Gauss theorem. 
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APPENDIX B: RECOVERING THE SHAPE 
EQUATION IN THE ABSENCE OF "ARTIFACT" 

The customary generic shape equation corresponding 
to our system without the "artifact" is the following co- 
variant equation (25j : 

^AH + H[H^ - K] = 0, (Bl) 

where A denotes the Laplacian, H the mean curvature, 
and K the extrinsic curvature. Once the Laplacian is 
expanded, the generic shape equation (IBip writes 

i(g^'5z - g""r%„„)a,if + H[H^-K]= 0, (B2) 

or, in isometric azimuthal coordinates {u,ip), 

\ e-^'^^'^'^duuHiu) + H{u) [h'^{u) - K{u)\ = 0. (B3) 

Substituting formulas (|A57P and (|A58P into previous 
Eq. dm leads to 

{pu — cosfi(u)) {duu^{u) — sinr2(u) cos51(m)) 

+ \{du^{u) - &mn{u}f {duVL{u) + smn{u)) = 0. 

(B4) 

On the other hand, our generic shape equation (PD|) 
without the "artifact" reads 

+ ^0 cos Q,{u) 
+ uq\^z{u) cos ^{u) + r{u) sin ^l{u)\ —0. (B5) 



In order to compare the shape equation (jB5[) with the 
general shape equation (|B4|) . let us first transform it in 
a form free of radius r{u), height z{u), and constants of 
integration fiQ and vq as follows. Since formulas (jA51a[) . 
(|A22|) . and (|A52|) yield the equality 

duz{u) cos^liu) + duriu) sii\il[u) = 0, (B6) 

the differentiation of Eq. (|B5p readily gives 



duu^{u) — sin r2(tt) cosr2(w) = — sinf2(w) 

+ t'o [— sin Q.{u) + r{u) cos . (B7) 

By applying the operator {pu — cosr2(it)) to both sides 
of previous equality (|B7|) . we obtain 

(pu — cos VL{u)) {puu^{u) — sin il{u) cos ~ 
— [no cos il{u) + vq [z{u) cos il{u) + r{u) sin ) 

- sinn(it)), (B8) 

after simplification of the right-hand side. Then the 
insertion of formula (|B5P into previous Eq. (jBSP gives 
Eq. (lB4l) . 

In summary, as announced, the two shape equations 
(jBSP and (jB4|) are exactly the same. 
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